Marginal Fermi liquid behavior from 2d Coulomb interaction 
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A full, nonperturbative renormalization group analysis of interacting electrons in a graphite layer 
is performed, in order to investigate the deviations from Fermi liquid theory that have been observed 
in the experimental measures of a linear quasiparticle decay rate in graphite. The electrons are 
coupled through Coulomb interactions, which remain unscreened due to the semimetallic character 
of the layer. We show that the model flows towards the noninteracting fixed-point for the whole 
range of couplings, with logarithmic corrections which signal the marginal character of the interaction 
—} ' separating Fermi liquid and non-Fermi liquid regimes. 
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During recent years there has been important progress in understanding the properties of quantum electron liquids 
I i in dimension D < 3. One of the most fruitful approaches in this respect springs from the use of renormalization group 
+3 ' (RG) methods, in which the different liquids are characterized by several fixed-points controlling the low-energy 
^ ■ properties. The Landau theory of the Fermi liquid in dimension D > 1 can be taken as a paradigm of the success 
of this program. It has been shown that, at least in the continuum limit, a system with isotropic Fermi surface and 

S regular interactions is susceptible of developing a fixed-point in which the interaction remains stable in the infrareda. 
The question of whether different critical points may arise at dimension D = 2 is now a subject of debateou. From 
the perspective of the RG approach, one of the premises leading to the Fermi liquid fixed-point should be relaxed in 
<**] ■ order to flow to a different universality class. In the case of models proposed to understand the electronic properties of 
copper-oxide superconductors, the high anisotropy of the Fermip*urfacd3 may play an important role in the anomalous 
behavior of the normal as well as of the superconducting state I. On the other hand, a possible source of non- Fermi 
liquid behavior may arise in systems with singular interactions!! In the case of the Coulomb interaction screened by 
t-H , the Fermi sea, a solution by means of bosonization methods has shown that no departures from Fermi liquid behavior 
arise at D = 2 and 3Q. It has been also shown for the conventional screened interaction that only potentials as singular 
as V(q) ~ l/\q\ 2D ~ 2 can lead to a different electron liquiaSlld. The system of electrons with magnetic interactions is 
quite different in that respect. It is known that in this case the system shows non-Fesrai liquid behavior for D < 3, 
| manifested in properties like the specific heat or anomalous electron field dimensionaljil3. 

In the present work, we address the applicability of the notion of Fermi liquid fixed-point to two-dimensional systems 
with unscreened Coulomb interaction. The absence of screening requires the vanishing of the density of states at the 
^\ , Fermi level. Semimetals show this behavior, while retaining a gapless electronic spectrum— The existence of a well 
defined continuum at low energies permits the existence of non trivial scaling propertieslij. The most remarkable 
example of this—kind is given by the two-dimensional sheet of graphite, which has a vanishing density of states at 
Q the Fermi levellj. Recent photoemissionj-experiments in graphite, at intermediate energies, show a decay rate of 
quasiparticles proportional to their energyt3. This represents a clear deviation with respect to the behavior in metals, 
which follow the conventional Fermi liquid picture with quasiparticle lifetimes propoportional to the inverse of the 
energy square, with, at most, logarithmic corrections. A description in terms of an effective field theory model has 
q , shown that the electronic interactions within the graphite layers are mainly responsible for the anomalous properties 
' measured in the experimenttZI. 
. ! We apply RG techniques to investigate whether the mentioned anomalous behavior can be understood as a marginal 
deviation from Fermi liquid theory, or rather it points towards a different universality class realized in the graphite 
sheet. We recall that the low-energy electronic excitations of the latter at half- filling are concentrated around two 
Fermi points at the corners of the hexagonal Brillouin Zone, where the dispersion relation is well approximated by 
two cones in contact at the apex. The effective field theory is given therefore by a pair of Dirac fermions, with a 
Coulomb potential that remains unscreened due to the vanishing density of states at the Fermi points. The effective 
hamiltonian can be written in the formllj 

H = -iv F ( d 2 rV + (r)cr ■ V*(r) + — / d 2 r x ( d 2 r 2 *+(ri)*(ri) 1 r * + (r 2 )*(r 2 ) (1) 

J 87T J J |rx-r 2 | 

where ^(r) is a two-dimensional Dirac spinor and cr = (a xi a y ). Such effective field theory provides a good starting 
point for a RG analysis since, given that the scaling dimension of the W(r) field is —1 (in length units), the four-fermion 
Coulomb interaction turns out to be scale invariant, at this level, with a dimensionless coupling constant e 2 . 
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In order to address the existence of a new universality class, besides the trivial noninteracting phase, a nonpertur- 
bative approach has to be adopted, since such new phase can only be revealed by a nontrivial fixed-point in coupling 
constant space. In the following, we will implement a GW approximation in the computation of the self-energy prop- 
erties. This is more easily achieved in the present model by replacing the four-fermion term in (Q) by the interaction 
with an auxiliary scalar field used to propagate the Coulomb interaction. The effective hamiltonian can be rewritten 
in the form 

H = -iv F J d 2 r^+{r)cr • W(r) +e J d 2 r#+(r)*(r) 0(r) (2) 

where the scalar field <fi(r) has the propagator 

i<T#r, t) 4,(r', f )) = ±5(t (3) 

In this framework we will introduce the GW approximation by taking into account the quantum corrections to the 
<f> propagator due to particle-hole excitations of the Fermi sea. This kind of approximation has proven to be adequate 
to the description of the crossover from Fermi liquid to Luttinger liquid behavior upon lowering the dimension from 
D = 2 to 1, capturing the relevant physical processes in the electron systemll3. Therefore, it seems also appropriate to 
uncover any possible fixed-point, different from that of Fermi liquid theory, in the case of the system with unscreened 
Coulomb interaction. 

The perturbative analysis of our model shows, in fact, the existence of a free fixed-point that is stable in the infrared 
limit. This can be understood from the nontrivial scaling of the model with respect to variations of the bandwidth 
cutoff E c , that is needed to regulate the divergent contribution of virtual processes to observable quantities. According 
to the RG point of view, a reduction of the cutoff E c implements a partial integration of the high-energy electron 
modes, that renormalize in this way the value of the effective parameters in the low-energy theory. The electron charge 
e is not renormalized in the present model, but the Fermi velocity vp is renormalized to first order in perturbation 
theory by a self-energy correction of the fornO 

e 2 

E(k,w fc ) « — <r • fclogi; c (4) 

Thus, in the perturbative regime the Fermi velocity Vp grows steadily as the cutoff E c is reduced, and the effective 
coupling constant e 2 /vp flows to zero in the low-energy effective theory. 

The most interesting point, however, concerns the analysis of the model away from the perturbative regime. With 
regard to the direct application to the graphite system the weak coupling results are of little use, since the bare 
coupling in the graphite sheet has an estimated value e 2 /vp ~ 10. Some possibly relevant effects like, for instance, 
the renormalization of the quasiparticle weight have to be consistently understood in a nonperturbative framework. 

As stated above, the polarization tensor does not show any divergence with respect to the bandwidth cutoff E c 
and, at the one-loop level, it is given by 

e 2 k 2 

jn(k,w fe ) = i- ; (5) 

The dressed propagator of the interaction in the RPA is given by 

— % 

{<t>(\L,Wk)<f>(-k,-U}k)) = „„ , , „ 2 k5 (6) 
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By using the dressed propagator in the computation of the self-energy one is able to perform a partial sum 
of perturbation theory, in which one takes into account the set of most singular diagrams with regard to the bare 
interaction ~ l/|k|. In this kind of GW approximation and for the model with conical dispersion relation, we have 

. , . 2 f d 2 p dw u; k - uj + v F (T-(k - p) -i 

12j(K, OJk) = *2e / 5— -5 — ; — 5 2 (7) 

J (2tt) 2 2tt v 2 F (k - p) 2 - K - c^) 2 2 |p| + f . 2 p2 2 ■ 

The imaginary part of the self-energy coming from (Q) has been computed elsewhereQ, and it has been shown to 
have a linear dependence on quasiparticle energy, consistent with the measured quasiparticle lifetimes in graphite. In 
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this paper we are interested in the computation of the real part of E(k, u;fc), which provides information about the 
nontrivial scaling of the quasiparticle weight and the Fermi velocity in the low-energy limit. 

The terms linear in u>k and k in the self-energy (Q) display a logarithmic dependence on the high-energy cutoff 
E c . This can be determined in the following way. Since we arc interested in the real part of E we may perform the 
analytic continuation lu — > i 57. We end up with the expression: 

e 2 f d 2 p duJ 1 uJk — SJ + cr-^vpk — p) 



e 2 f 
vf J 



(27r) 2 27r |p| {v F k - p) 2 + (ZJ k - ZJV 



. e 2 e 2 f d 2 p duo Vp 2 +u 2 ~ jt^\p\ uJ k - U + cr-(v F k - p) 
'v f 16vfJ (2ir) 2 2it ap 2 +H 2 (u^k — p) 2 + (o7 fc — ul) 2 



where a = 1 — (e 2 /(16vi?)) 2 - A singularity at a = appears in the above expression, whose role has to be clarified 
since there is no sign of a particular feature at e 2 /(16vp) — 1 in the original expression. We will perform the above 
integrals taking the value of a > 0, but it will become clear at the end that the results can be continued smoothly to 
the strong coupling regime e 2 /(16vf) > 1. 

Upon integration of the frequency from — oo to +oo and placing the bandwidth cutoff E c in momentum space, 
Vi?|p| < E c , the coefficients of the logarithmically divergent contributions can be computed in terms of elementary 
functions of g = e 2 / (16vf). The renormalization of the electron propagator turns out to be given by 

G = G^ _S 

= Z^ 1 (ujk — Vpu-k) 

-Z-^" (g 2 + (2 - g 2 ) ( 1 - ^=£^ ) lo g E c - ( l=££ - 1^ log^ 



g^l-g 2 / Jr g V y/i - g? 



+Z- 1 v F er-kX ( 1 - g2 arcsin 5 | log^ c - Z^vpa-k-- (l - J\ - g 2 ) \ogE c (9) 

7T \ 9 J 7I"5 V ' 

where Z 1 ' 2 represents the scale of the bare electron field compared to that of the cutoff-independent electron field 

VbareiEc) = Z^H (10) 

In the RG approach, wc require the cutoff-independence of the renormalized Green function, since this object leads 
to observable quantities in the quantum theory. For this purpose, Z and vp have to be understood as cutoff-dependent 
effective parameters, that reflect the behavior of the quantum theory as E c — > and more states are integrated out 
from high-energy shells of the band. We get the RG flow equations 

„ d , „ . 8 / 2 — q 2 arccos q \ 8 1 , . 

E c —)ogZ{E c ) = — i \2+ 7= = +— 11 

dE c it 2 y g y/\ - g 2j irg 

t-, d ( n \ 8 / arccos g \ 4 1 . „. 

E c —v F (E c ) = -vf 1 + —7== + -vf- 12 

dE c IT 2 y g^Jl -glj 7T g 

Given that the electron charge e is not renormalized, we may write down the flow equation for the effective coupling 
constant g — e 2 /(16i>F) 

d , . 8 1 arccos q \ 4 . , 

E c —g(E c ) = -j [g + 1 - - 13 

dE c n 2 y y/i - g 2 J 7T 

In the weak coupling regime, one may check that the renormalization of both the Fermi velocity and the electron 
wavefunction takes place in the expected direction. The quasiparticle weight Z at small g is smaller than the bare 
value measured before integration of high-energy modes. The Fermi velocity vp flows to higher values in the infrared, 
and the density of states around the Fermi energy decreases, as a consequence of screening effects. This ensures the 
consistency of the weak coupling phase, where the results of perturbation theory become increasingly reliable in the 
low-energy limit. However, the most important point concerns the possible existence of a different phase at large 
values of g. In this respect, the flow equations (jll) and ( |l3| ) can be analytically continued to values g > 1, by simple 
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use of the formula arccosz = ilog(z — iyl — z 2 ). The flows of the coupling constant and the electron wavefunction 
are then differentiable across 5 = 1, which shows that the apparent singularity at this point has no physical meaning. 

The right-hand-side of Eq. (|l3|) is a monotonous function of g, taking into account the mentioned analytic continu- 
ation. This means that there is no phase different from that of the perturbative regime, and that the strong coupling 
regime is connected to it through RG transformations. The RG flow represented in Fig. [I] shows that the perturbative 
regime is attained at low-energies, starting from fairly large bare values of the coupling constant. 

The present analysis is relevant to the phenomenology of the graphite layers. It shows that, even in such a 
system with unconventional quasiparticle lifetimes, the low-energy behavior is governed by a fixed-point which can be 
described as a Fermi liquid, as Z tends to a constant in the infrared, unlike for the line of nontrivial fixed-points which 
characterize Luttinger liquids. The low-energy scaling of the quasiparticle weight is represented in Fig. |^. Hence 
we can assert that, although the quasiparticles decay according to the experiment with a rate proportional to their 
energy, rather than to their energy square, the notion of well-defined quasiparticles with a finite weight at the Fermi 
level still holds in the system. 

Our study stresses the anomalous screening properties of the Coulomb interaction in low-dimensional systems. This 
fact has been also put forward recently in a different framework, pointing out that in one and two-dimensional systems 
the screening of the long-range interactions goes in the direction of reducing the electron correlational. In the RG 
approach we see how this effect arises under the form of a renormalization of vp. It is remarkable that such nontrivial 
scaling of the Fermi velocity in the infrared operates in the two-dimensional as well as in the one-aLLm£nsional system 
with Coulomb interaction! 2 .!], and it seems to be also present in systems with magnetic interactional^. 

The analysis presented here is relevant to the controversy about the existence of systems with Fermi liquid behavior 
in two dimensions. Our approach is completely rigorous in the limit of a very large number N of electron flavors, 
where the bubble summation implicit in the RPA becomes exact. The leading order in the 1/N expansion provides, 
on the other hand, the fixed-point and the exact anomalous dimension of the electron field in the Luttinger modelr 9 !! 21 !. 
Therefore, this approximation should suffice to describe such kind of nontrivial fixed-point, in case that it were present 
in the model. 

To summarize, semimetals described by the two-dimensional Dirac equation, such as a graphite layer, show signif- 
icant differences with respect to the properties of standard Fermi liquids with (screened) Coulomb interactions. In 
the present problem, the quasiparticle lifetime goes like ~ w" 1 , while the enhancement of the Fermi velocity implies 
the vanishing of the effective coupling in the infrared. We expect the same behavior in three-dimensional zero-gap 
semiconductors, which are also described by an effective Dirac equation. 

In the context of more general long-range interactions, the system with a Dirac sea behaves also differently with 
respect to the conventional Fermi sea, as in the former an interaction V(q) ~ l/|(7| 1+e already departs from the Fermi 
liquid universality class for e > 0. In our RG framework the interaction becomes relevant no matter how small e may 
be, and a nontrivial fixed-point can be found away from the origin within the e expansion. 

The logarithmic behavior that we have found affects also some thermodynamic quantities like the specific heat or 
the susceptibility, which pick up logarithmic corrections as T — ► 0. These are the signature of the marginal character 
of the interaction, which has in our model the precise degree of singularity to separate regimes with Fermi liquid 
(e < 0) and non- Fermi liquid behavior (e > 0). 
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FIG. 1. Flow of the coupling constant for different bare values. 



FIG. 2. Wavefunction renormalization for bare coupling constant g — 5 (thick line) and g = 1 (thin line). 
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